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We study the influence of laser radiation on the electron transport through a molecular wire 
weakly coupled to two leads. In the absence of a generalized parity symmetry, the molecule rectifies 
the laser induced current resulting in directed electron transport without any applied voltage. We 
consider two generic ways of dynamical symmetry breaking: mixing of different harmonics of the 
laser field and molecules consisting of asymmetric groups. For the evaluation of the nonlinear current, 
a numerically efficient formalism is derived which is based upon the Floquet solutions of the driven 
molecule. This permits a treatment in the non-adiabatic regime and beyond linear response. 
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I. INTRODUCTION 

During the last years, we experienced a wealth 
of experinaetital activity in the field of molecular 
electronics Its technological prospects for nano- 
circuitstl have created broad interest in the conductance 
of moleculea^tached to metal surfaces or tips. In recent 
experimentaj'tlErLl weak tunneling currents through only a 
few or even single molecules coupled by chemisorbed thiol 
groups to the gold surface of leads has been achieved. 
The experimental development is accompanied by an in- 
creasing theoretical interest in the transport properties 
of such systems.l3'3 An intriguing challenge presents the 
possibility to control the tunneling current through the 
molecule. Typical energy scales in molecules are in the 
optical and the infrared regime, where today's laser tech- 
nology provides a wealth of coherent light sources. Hence, 
lasers represent an inherent possibility to control atoms 
or molecules and to direct currents through them. 

A widely studied phenomenon in extended, strongly 
driven ..dijivpn —systems is the so-termed ratchet 
effect ,E£lllij'Ejll30 originally discovered and investigated 
for overdamped classical Brownian motion in periodic 
nonequilibrium systems in the absence of reflection sym- 
metry. Counterintuitively to the second law of thermody- 
namics, one then observes a directed transport although 
all acting forces possess no net bias. This effect has been 
established as well within the regime-of dissipative, in- 
coherent quantum Brownian motionJl3 A related effect 
is found in the overdamped limit of dissipative tunnel- 
ing in tight-binding lattices. Here the spatial symme- 
try is typically preserved and the nonvanishing transport 
is brought about by harmonic mixing, of a driving field 
that includes higher harmonics. t£lEZlE3 For overdamped 
Brownian motion, both phenomena can be understood, 
in terms of breaking a generalized reflection symmetry.t2l 

Recent theoretical descriptions of m olecu lar conductiv- 
ity are based on a scattering approachB3E]l Alternatively, 
one can assume that the underlying transport mechanism 
is an electron transfer reaction and that the conductivity 



can be derived from the corresponding reaction rate.u 
This analogy leads to a connection between electron 
transfer rates in a donor-acceptor system and conduction 
in the same system when-operating as a molecular wire 
between two metal leads.t3 Within the high-temperature 
limit, the electron transport x)EL.ttLft_syire can be described 
by inelastic hopping events E'CJ'oES For a more quanti- 
tative ah initio analysis, the molecular orbitals may be 
taken from electronic structure calculations .l3 

Isolated atoms and molecules in strong oscillat- 
ing fiplHg hp'^'e i^^'^f"., •'Y''^'^b^ studied within a Floquet 
formalisinEZl'tjO'L2lEirE3 and many corresponding theo- 
retical techniques have been developped in that area. 
This suggests the procedure followed in Rcf. |3^: Making 
use of these Floquet tools, a formalism for the transport 
through time-dependent quantum systems has been de- 
rived that combines Floquet theory for a driven molecule 
with the many-particle description of transport through 
a system that is coupled to ideal leads. This approach 
is devised, much in the spirit of the Floquet-Markov 
theoryoEj for driven dissipative quantum systems. It as- 
sumes that the molecular orbitals that are relevant for the 
transport are weakly coupled to the contacts, so that the 
transport characteristics are dominated by the molecule 
itself. Yet, this treatment goes beyond the usual rotating- 
wave approximation as frequently employed, such as e.g. 
in Refs. ^5[|3^ . 

A time-dependent perturbative approach to the prob- 
lem of driven molepalaj' wires has recently been described 
by Tikhonov et aZ.E3't3 However, their one-electron treat- 
ment of this essentially inelastic transmission process 
cannot handle consistently the electronic populations on 
the leads. Moreover, while their general formulation is 
not bound to their independent channel approximation, 
their actual application of this approximation is limited 
to the small light-molecule interaction regime. 

With this work we investigate the possibilities for 
molecular quantum wires to act as coherent quantum 
ratchets, i.e. as quantum rectifiers for the laser-induced 
electrical current. In doing so, we provide a full account 
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of the derivation published in letter format in Ref. In 
Sec. II we present a more detailed derivation of the Flo- 
quet approach to the transport through a periodically 



[II to in- 



driven wire. This formalism is employed in Sec. 
vestigate the rectification properties of driven molecules. 
Two generic cases are discussed, namely mixing of differ- 
ent harmonics of the laser field in symmetric molecules 
and harmonically driven asymmetric molecules. We fo- 
cus thereby on how the symmetries of the model sys- 
tem manifest themselves in the expressions for the time- 
averaged current. The general symmetry considerations 
of a quantum system under the influence of a laser field 
are deferred to the Appendix ^ 
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FIG. 1: Level structure of a molecular wire with A'' = 8 
atomic sites which are attached to two leads. 



II. FLOQUET APPROACH TO THE 
ELECTRON TRANSPORT 

The entire system of the driven wire, the leads, and the 
molecule-lead coupling as sketched in Fig. |l| is described 
by the Hamiltonian 



H{t)^H„i,,{t)+Hi 



leads 



wirc-lcads ■ 



(1) 



The wire is modeled by N atomic orbitals |n), n = 
1, . . . , A'^, which are in a tight-binding description cou- 
pled by hopping matrix elements. Then, the correspond- 
ing Hamiltonian for the electrons on the wire reads in a 
second quantized form 



(2) 



where the fermionic operators c„, c]j annihilate, respec- 
tively create, an electron in the atomic orbital \n) and 
obey the anti-commutation relation [c„,Cjj,]+ — 6n.n'- 
The influence of the laser field is given by a periodic time- 
dependence of the on-site energies yielding a single parti- 
cle Hamiltonian of the structure i?„„/(t) = Hnn'{t + T), 
where T = 2Tr/fl is determined by the frequency ft of the 
laser field. 

The orbitals at the left and the right end of the 
molecule, that we shall term donor and acceptor, |1) and 
I TV), respectively, are coupled to ideal leads (cf. Fig. |^) 
by the tunneling Hamiltonians 



;-lcads 



VqR cIj^Cn) 



H.c. 



(3) 



The operator CqL {cqr) annihilates an electron on the 
left (right) lead in state Lq (Rq) orthogonal to all wire 
states. Later, we shall treat the tunneling Hamiltonian 
as a perturbation, while taking into account exactly the 
dynamics of the leads and the wire, including the driving. 

The leads are modeled as non-interacting electrons 
with the Hamiltonian 



leads 



qL ^'qL<-qL "t" tqR <^qR<^qR 



(4) 



A typical metal screens electric fields that have a fre- 
quency below the so-called plasma frequency. Therefore, 
any electromagnetic radiation from the optical or the in- 
frared spectral range is almost perfectly reflected at the 
surface and will not change the bulk properties of the gold 
contacts. This justifies the assumption that the leads are 
in a state close to equilibrium and, thus, can be described 
by a grand-canonical ensemble of electrons, i.e. by a den- 
sity matrix 

fflcads.cq OC CXp [-(i/lcads - /-JlA^L - AiflA^i?^/^^^] , (5) 

where ^l/r. the electro-chemical potentials and 



L/R. 



q ^Il/r'^QL/R 



the electron numbers in the 



left/right lead. As a consequence, the only non-trivial 
expectation values of lead operators read 



(6) 



where e^L is the single particle energy of the state qL 
and correspondingly for the right lead. Here, f{x) = 
(1 -I- e^/ksT^j-i (jguotes the Fermi function. 



A. Time-dependent electrical current 

The net (incoming minus outgoing) current through 
the left contact is given by the negative time derivative 
of the electron number in the left lead, multiplied by the 
electron charge — e, i.e. 

n IP 

I^{t) = e-{N^),^-{[H{t),NL]),. (7) 

Here, the angular brackets denote expectation values at 
time i, i.e. {0)t = Tr[Op(<)]. The dynamics of the den- 
sity matrix is governed by the Liouville-von Neumann 
equation ihg{t) = [H{t), g{t)] together with the factoriz- 
ing initial condition gito) = gwirdto) (8) fteads^eq- For the 
Hamiltonian (|^), the commutator in Eq. (^ is readily 
evaluated to 



hit) 



qL{CqLCi)t 



(8) 
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To proceed, it is convenient to switch to the interaction 
picture with respect to the uncoupled dynamics, where 
the Liouville-von Neumann equation reads 

ih-^g{t, to) = [ffwiro-lcads(i, ^o), ^(i, to)]. (9) 

The tilde denotes the corresponding interaction picture 
operators, X{t,t') = U^{t,t') X{t)Uoit,t'), where the 
propagator of the wire and the lead in the absence of the 
lead-wire coupling is given by the time-ordered product 

Uo{t,t') = Texp J dt" [H^ireit") + i/lcads]) • 

(10) 

Equation is equivalent to the integral equation 

g(t,to) ^ gito,to) - T / dt'[H^iro-\eiids{,t' ,to), g{t' ,to)]. 

(11) 

Inserting this relation into Eq. (^, we obtain an ex- 
pression for the current that depends on the density of 
states in the leads times their coupling strength to the 
connected sites. At this stage it is convenient to intro- 
duce the spectral density of the lead-wire coupling 

rL/R{e) = Y E Kl/r\'K^ ~ ^.l/r), (12) 

which fully describes the leads' influence. If the lead 
states are dense, r^/^(e) becomes a continuous function. 
Since we restrict ourselves to the regime of a weak wire- 
lead coupling, we can furthermore assume that expecta- 
tion values of lead operators are at all times given by their 
equilibrium values (^. Then we flnd after some algebra 
for the stationary (i.e. for to — * — oo), time- dependent net 
electrical current through the left contact the result 

oo 

hit) ^^ReJdTj deVUe) e''^/^{{c\ 5i(t, t - t))^_^ 



- [cl.ci(i,i - T)]+f{e - ^l)}- 

(13) 

A corresponding relation holds true for the current 
through the contact on the right-hand side. Note that the 
anti-commutator [c|,ci(t,t — t)]+ is in fact a c-number 
(see below). Like the expectation value (c|ci(t, i — 
t))j_^, it depends on the dynamics of the isolated wire 
and is influenced by the external driving. 

It is frequently assumed that the attached leads can be 
described by a one-dimensional tight-binding lattice with 
hopping matrix elements A'. Then, the spectral densi- 
ties Ti^/ji{c) of the lead-Jiz.ire couphngs are given by the 
Anderson-Newns model,E3 i.e. they assume an elliptical 
shape with a band width 2A'. However, because we are 



mainly interested in the behavior of the molecule and not 
in the details of the lead-wire coupling, we assume that 
the conduction band width of the leads is much larger 
than all remaining relevant energy scales. Consequently, 
we approximate in the so-called wide-band limit the func- 
tions T i^/ji{e) by the constant values F^/^. The first con- 
tribution of the e-integral in Eq. (|l^) is then readily eval- 
uated to yield an expression proportional to 5{t) . Finally, 
this term becomes local in time and reads eFi(c|ci)^. 

B. Floquet decomposition 

Let us next focus on the single-particle dynamics of 
the driven molecule decoupled from the leads. Since its 
Hamiltonian is periodic in time, i?„„' (t) = Hnn' (t + 
T), we can solve the corresponding time-dependent 
Schrodinger equation within a Floquet approach. This 
means that we make use of the fact-^hat-thereL exists a 
complete set of solutions of the forn£3'BBEEa 

|*a(t)) =e-'^^°*/''|$„(t)), |$„(i)) = |$„(t + r)) (14) 

with the quasienergies £„. Since the so-called Floquet 
modes |^a(i)) obey the time-periodicity of the driving 
field, they can be decomposed into the Fourier series 

|1'a(t)) =Ef^"'''"l*"-fc)- (15) 
fc 

This suggests that the quasienergies come in classes, 

Ca.k^ea + knn, /c = 0,±l,±2,... , (16) 

of which all members represent the same solution of 
the Schrodinger equation. Therefore, the quasienergy 
spectrum can be reduced to a single "Brillouin zone" 
—Ml/2 < e < Ml/2. In turn, all physical quantities 
that are computed within a Floquet formalism are inde- 
pendent of the choice of a specific class member. Thus, 
a consistent description must obey the so-called class in- 
variance, i.e. it must be invariant under the substitution 
of one or several Floquet states by equivalent ones, 

ea, |$«(i)) ^ e« + k^rin, e'*^""*|$,(i)), (17) 

where fci, . . . , fcjv are integers. In the Fourier decomposi- 
tion (p^, the prefactor exp(ifcQJ7t) corresponds to a shift 
of the side band index so that the class invariance can be 
expressed equivalently as 

Ca, > + kaTlQ., |<I>„,fc+fc„}. (18) 

Floquet states and qua|Sifjn er m pfica n. be obtained from 
the quasienergy equationEZrEflatj'EiLJ 

(El^^)^"»'(*)("'|-*?i^)l'i'a(i)) =ea|'i'a(t))- (19) 

n,n' 

A wealth of methods for the solution of Ais, eigenvalue 
problem can be found in the literature. EirE3 One such 
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method is given by the direct numerical diagonaUzation 
of the operator on left-hand side of Eq. (|l^) . To account 
for the periodic time-dependence of the \^a{t)), one has 
to extend the original Hilbert space by a T-periodic 
time coordinate. For a harmonic driving, the eigenvalue 
problem (^9|) is band-diagonal and selected eigenvalues 
and eigenvectorsLxiaiL be computed by a matrix-continued 
fraction schemeOo 

In cases where many Fourier coefficients (in the present 
context frequently called "sidebands" ) must be taken into 
account for the decomposition (^5|), direct diagonaUza- 
tion is often not very efficient and one has to apply more 
elaborated schemes. For example, in the case of a large 
driving amplitude, one can treat the static part of the 
Hamiltonian as a perturbation.Eaoc3 The Floquet states 
of the oscillating part of the Hamiltonian then form an 
adapted basis set for a subsequently more efficient nu- 
merical diagonalization. 

A completely different strategy to obtain the Floquet 
states is to propagate the Schrodinger equation for a 
complete set of initial conditions over one driving pe- 
riod to yield the one-period propagator. Its eigenvalues 
represent the Floquet states at time t = 0, i.e., |$q(0)). 
Fourier transformation of their time-evolution results in 
the desired sidebands. Yet another, very efficieiirf: propa- 
gation scheme is the so-called (t, t')-formalism.Cj 

As the equivalent of the one-particle Floquet states 
\^a{t)), we define a Floquet picture for the fermionic 
creation and annihilation operators cjj, c„, by the time- 
dependent transformation 



Cait) = ^($„(t)|n) c„. 



The inverse transformation 



= ^(n|$„(t)}c„(i) 



(20) 



(21) 



follows from the mutual orthogonality ancLtke complete- 
ness of the Floquet states at equal timesH'o Note that 
the right-hand side of Eq. ( pT| ) becomes i-independent 
after the summation. In the interaction picture, the op- 
erator Ca{t) obeys 



dait,t') = ul{t,t')c^{t)Ua{t,t') 



(22) 



This is easily verified by differentiating the definition in 
the first line with respect to t and using that \^a{t)) 
is a solution of the eigenvalue equation (|lj). The fact 
that the initial condition Ca{t',t') — Ca{t') is fulfilled, 
completes the proof. Using Eqs. ( ^ ) and (|2^), we are 
able to express the anti-commutator of wire operators at 



different times by Floquet states and quasienergies: 
[c„,,ct(t,t')]+ = $]e""(*-*')/''(n'l$„(i'))(<fa(Oki)- 

a 

(23) 

This relation together with the spectral decomposi- 
tion (|l^) of the Floquet states allows to carry out the 
time and energy integrals in the expression ( p^ for the 
net current entering the wire from the left lead. Thus, 
we obtain 



(24) 



with the Fourier components 



{'^a,k'+k"\l){M'^f3,k+k")Ral3.k' 

aPk'k" 
a.k' 

+ {<Pa,k'-k\l){ma,k'))f{ 

Here, we have introduced the expectation values 



E 



^—ikflt I) 



(25) 

(26) 
(27) 



The Fourier decomposition in the last line is possible be- 
cause all Rap{t) are expectation values of a linear, dis- 
sipative, periodically driven system and therefore share 
in the long-time limit the time-periodicity of the driving 
field. In the subspace of a single electron, Rap reduces 
to the density matrix in the basis of the Floquet states 
which has been used to describe dissipative driven quan- 
tum systems in Refs. |32l||,PPr 



C. Master equation 

The last step towards the stationary current is to find 
the Fourier coefficients i?a/3,fe at asymptotic times. To 
this end, we derive an equation of motion for the reduced 
density operator Pwire(i) = Tricads g(i) by reinserting 
Eq. (|ri|) into the Liouville-von Neumann equation (||). 
We use that to zeroth order in the molecule-lead cou- 
pling the interaction-picture density operator does not 
change with time, g(t — t, to) « g{t,to). A transforma- 
tion back to the Schrodinger picture results after tracing 
out the leads' degrees of freedom in the master equation 
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Qwii-c{t) — —^[Hmrc(t), gl„iro(0] ^ j Tlioads [-f^wiro-lcads, [^^wiro - leads — T,t), glwirc(i) ® £'lcads,cq]] • (28) 



Since we only consider asymptotic times to — > — oo, we 
have set the upper Umit in the integral to infinity. From 
Eq. ( p8|) follows directly an equation of motion for the 
Rap{t). Since all the coefficients of this equation, as well 
as its asymptotic solution, are T-periodic, we can split it 
into its Fourier components. Finally, we obtain for the 
Raf3,k the inhomogeneous set of equations 



(29) 



^ E ( E {^P,k'+k"\i)mp',k+k")Ra0'.k' 

k' 13' k" 

+ ^{^c.'.k'+k"\i)m a,k-\-k ")Ra'f3.k' 



a'k" 



- ($/3,fc'-fc|l)(l|$a,/c')/(ea,/c' - Mi) 

- ($/3^fc'|l)(l|*a,fe' + fc>/(e/3,fe' - ml) 

h same terms with the replacement 



For a consistent Floquet description, the current formula 
together with the master equation must obey class in- 
variance. Indeed, the simultaneous transformation with 
( |l8| ) of both the master equation ( p9| ) and the current 
formula ( p5| ) amounts to a mere shift of summation in- 
dices and, thus, leaves the current as a physical quantity 
unchanged. 

For the typical parameter values used below, a large 
number of sidebands contributes significantly to the 
Fourier decomposition of the Floquet modes |$a(t)). Nu- 
merical convergence for the solution of the master equa- 
tion (|2^), however, is already obtained by just using a few 
sidebands for the decomposition of Ra/3{t). This keeps 
the numerical effort relatively small and justifies a poste- 
riori the use of the Floquet representation (^l|). Yet we 
are able to treat the problem beyond a rotating-wave- 
approximation. 



through the left contact it is given by 

h ^ il ^ eTLj2[j2{^a.k'+k\i)mp.k')Rc.p,k 

ak I3k' 

- («'a,fe|l)(l|«'a,fe)/(ea,fc " Mi) 

(30) 

Mutatis mutandis we obtain for the time-averaged net 
current that enters through the right contact 

.k' + k\N){N\^p,k')Ral3.k 
ak Pk' (31) 

- ($a,fc|iV)(7V|$„,fe)/(e„,fe-Ai;^j . 

Total charge conservation of the original wire-lead 
Hamiltonian (^ of course requires that the charge on 
the wire can only change by current flow, amounting to 
the continuity equation Qwiro(t) — Ihii) + Init)- Since 
asymptotically, the charge on the wire obeys at most 
the periodic time-dependence of the driving field, the 
time-average of Qwiieit) must vanish in the long-time 
limit. From the continuity equation one then finds that 
II + Ir = 0, and we can introduce the time-averaged 
current 



I = Il = ~Ir. 



(32) 



For consistency, the last equation must also follow from 
our expressions for the average current ( |30| ) and ( ^ ) . In 
fact, this can be shown by identifying + Ir as the sum 
over the right-hand sides of the master equation ( p9| ) for 
a — P and /c = 0, 



Il + Ir = E 

a 

= 0, 

which vanishes as expected 



^(fa - e/3 + khn)Ra/3^k 



a=/3,k=0 



(33) 



D. Average current 

Equation ( p4|) implies that the current II {t) obeys the 
time-periodicity of the driving field. Since we consider 
here excitations by a laser field, the corresponding fre- 
quency lies in the optical or infrared spectral range. In 
an experiment one will thus only be able to measure the 
time-average of the current. For the net current entering 



E. Rotating-wave approximation 

Although we can now in principle compute time- 
dependent currents beyond a rotating-wave approxima- 
tion (RWA) , it is instructive to see under what conditions 
one may employ this approximation and how it follows 
from the master equation (^9|) . We note that from a com- 
putational viewpoint there is no need to employ a RWA 
since within the present approach the numerically costly 
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part is the computation of the Floquet states rather than 
the solution of the master equation. Nevertheless, our 
motivation is that a RWA allows for an analytical solution 
of the master equation to lowest order in the lead-wire 
coupling r. We will use this solution below to discuss 
the influence of symmetries on the F-dependence of the 
average current. 

The master equation ( p9| ) can be solved approximately 
by assuming that the coherent oscillations of all Ra(){t) 
are much faster than their decay. Then it is useful to 
factorize Rapit) into a rapidly oscillating part that takes 
the coherent dynamics into account and a slowly de- 
caying prefactor. For the latter, one can derive a new 
master equation with oscillating coefficients. Under the 
assumption that the coherent and the dissipative time- 
scales are well separated, it is possible to replace the 
time-dependent coefficients by their time-average. The 
remaining master equation is generally of a simpler form 
than the original one. Because we work here already with 
a spectral decomposition of the master equation, we give 
the equivalent line of argumentation for the Fourier co- 
efficients Ral3,k- 

It is clear from the master equation (E9|) that if 



(34) 



then the corresponding Rap^k emerge to be small and, 
thus, may be neglected. Under the assumption that the 
wire-lead couplings are weak and that the Floquet spec- 
trum has no degeneracies, the RWA condition (34) is well 
satisfied except for 



a — 13, k ~ 0, 



(35) 



i.e. when the prefactor of the l.h.s. of Eq. ( p4[ ) vanishes 
exactly. This motivates the ansatz 



Raf3,k — Pa ^a,f3 ^, 



Ofc,0, 



(36) 



which means physically that the stationary state consists 
of an incoherent population of the Floquet modes. The 
occupation probabilities Pa are found by inserting the 
ansatz (^^ into the master equation (^^ and read 



Eki^ 



) 



Thus, the populations are determined by an average over 
the Fermi functions, where the weights 



TL\ma,k)\' 



<k=^R\{ma,k)\' 



"a,k 
..N 



(38) 
(39) 



are given by the effective coupling strengths of the k-th 
Floquet sideband \^a.k) to the corresponding lead. The 
average current (|3^) is within RWA readily evaluated to 
read 



wX I.W, 



'RWA 



Sk'^k"iK, 



N 



■ W 



N 



(40) 



[/( 



) - /(Ca.fe - Mi)]- 



III. RECTIFICATION OF THE 
DRIVING-INDUCED CURRENT 

In the absence of an applied voltage, i.e. = fJ-R, 
the average force on the electrons on the wire vanishes. 
Nevertheless, it may occur that the molecule rectifies 
the laser-induced oscillating electron motion and conse- 
quently a non-zero dc current through the wire is estab- 
lished. In this section we investigate such ratchet cur- 
rents in molecular wires. 

As a working model we consider a molecule consisting 
of a donor and an acceptor site and N — 2 sites in be- 
tween (cf. Fig. 1^). Each of the sites is coupled to its 
nearest neighbors by a hopping matrix elements A. The 
laser field renders each level oscillating in time with a 
position dependent amplitude. The corresponding time- 
dependent wire Hamiltonian reads 



Hnn'{t) — — ^{S7un' + 1 + Sn+l,n') 

+ [En - a{t) 



(41) 



where Xn = (A^ + 1 — 2n)/2 is the scaled position of site 
the energy a[t) equals the electron charge multiplied 
by the time-dependent electrical field of the laser and the 
distance between two neighboring sites. The energies of 
the donor and the acceptor orbitals are assumed to be 
at the level of the chemical potentials of the attached 
leads, El — E^ ~ /i^ = ^r. The bridge levels i?„, 
n = 2, . . . , A — 1, lie Eb above the chemical potential, as 
sketched in Fig. |l|. Later, we will also study the modified 
bridge sketched in Fig. ^ below. We remark that for 
the sake of simplicity, intra-atomic dipole excitations are 
neglected within our model Hamiltonian. 

In all numerical studies, we will use the hopping matrix 
element A as the energy unit; in a realistic wire molecule, 
A is of the order 0.1 eV. Thus, our chosen wire-lead hop- 
ping rate F = 0.1A/?i yields eF — 2.56 x 10"^ Ampere 
and f2 = 3A/?i corresponds to a laser frequency in the in- 
frared. Note that for a typical distance of SA between two 
neighboring sites, a driving amplitude A = A is equiva- 
lent to an electrical field strength of 2 x 10^ V/cm. 



A. Symmetry 

It is known frornJ|he study of deterministically rocked 
periodic potemtialgcZl and of overdamped classical Brow- 
nian motiorJla that the symmetry of the equations of mo- 
tion may rule out any non-zero average current at asymp- 
totic times. Thus, before starting to compute ratchet cur- 
rents, let us first analyze what kind of symmetries may 
prevent the sought-after effect. Apart from the principle 
interest, such situations with vanishing average current 
are also of computational relevance since they allow to 
test numerical implementations. 

The current formula (25) and the master equation (B9) 



contain, besides Fermi factors, the overlap of the Flo- 
quet states with the donor and the acceptor orbitals |1) 
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and \N). Therefore, we focus on symmetries that re- 
late these two. If we choose the origin of the position 
space at the center of the wire, it is the parity transfor- 
mation V : X —X that exchanges the donor with the 
acceptor, |1) |A^). Since we deal here with Floquet 
states I respectively with their Fourier coefficients 

|$a,fc), we must take also into account the time t. This 
allows for a variety of generalizations of the parity that 
differ by the accompanying transformation of the time 
coordinate. For a Hamiltonian of the structure (41), 
two symmetries come to mind: a{t) = —a{t + ir/n) 
and a{t) = —a{—t). Both are present in the case of a 
purely harmonic driving, i.e. a{t) cx sin(f2t). We shall 
derive their consequences for the Floquet states in the 
Appendix ^ and shall only argue here why they yield a 
vanishing average current within the present perturbative 
approach. 



1. Generalized parity 

As a first case, we investigate a driving field that obeys 
a(t) = —a{t + Tr/rt). Then, the wire Hamiltonian (^ ) is 
invariant under the so-called generalized parity transfor- 
mation 



5, 
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{x,t) ^ {^x,t + n/n). 



(42) 



Consequently, the Floquet states are either even or odd 
under this transformation, i.e. they fulfill the relation 
(A5), which reduces in the tight-binding limit to 



(43) 



where a a = ±1, according the generalized parity of the 
Floquet state |$q(<)). 

The average current / is defined in Eq. (^2|) by the 
current formulae (^o|) and (|l|) together with the mas- 
ter equation (|2^). We apply now the symmetry relation 
( ^ ) to them in order to interchange donor state |1) and 
acceptor state |A^). In addition we substitute in both 
the master equation and the current formulae Ra0,k by 

Raf3,k — <^a<^p{~'^)'^ Raf},k- The rcsult is that the new ex- 
pressions for the current, including the master equation, 
are identical to the original ones except for the fact that 
Ilt^l and In^Tfi are now interchanged (recall that we 
consider the case = Mi?,)- Therefore, we can conclude 
that 



Ir 



(44) 



which yields together with the continuity relation ( |33| ) a 
vanishing average current 7 = 0. 



2. Time-reversal parity 

A further symmetry is present if the driving is an odd 
function of time, a(t) = —a{—t). Then, as detailed in 



the Appendix]^, the Floquet eigenvalue equation (19) i 
invariant under the time-reversal parity 
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(a>,a:,i)^($*,-x,-i). 



(45) 



i.e. the usual parity together with by time-reversal and 
complex conjugation of the Floquet states The con- 
sequence for the Floquet states is the symmetry relation 
(|A7| ) which reads for a tight-binding system 



(46) 



Inserting this into the current formulae (^0|) and ( |3l| ) 
would yield, if Rap,k were real, again the balance condi- 
tion (|44| ) and, thus, a vanishing average current. How- 
ever, the Ra/3,k are in general only real for F^ = F^ = 0, 
i.e. for very weak coupling such that the condition (|3^ ) 
for the applicability of the rotating-wave approximation 
holds. Then, the solution of the master equation is dom- 
inated by the RWA solution (|3^), which is real. In the 
general case, the solution of the master equation ( p9) i s 
however complex and consequently the symmetry (4q) 
does not inhibit a ratchet effect. Still we can conclude 
from the fact that within the RWA the average current 
vanishes, that / is of the order F^ for F — > 0, while it is 
of the order F for broken time-reversal symmetry. 



B. Rectification from harmonic mixing 



The symmetry analysis in Sec. [HA explains that a 
symmetric bridge like the one sketched in Fig. ^ will not 
result in a average current if the driving is purely har- 
monic since a non-zero value is forbidden by the general- 
ized parity (^2|) . A simple way to break the time-reversal 
part of this symmetry is to add a second harmonic to the 
driving field, i.e., a contribution with twice the funda- 
mental frequency O, such that it is of the form 



a{t) = Ai sm{nt) + Ai sin(2ilt 



(47) 



as sketched in Fig. |^. While now shifting the time t by 
a half period Tr/fi changes the sign of the fundamental 
frequency contribution, the second harmonic is left un- 
changed. The generalized parity is therefore broken and 
we find generally a non- vanishing average current. 

The phase shift plays here a subtle role. For = (or 
equivalently any multiple of tt) the time-reversal parity 
is still present. Thus, according to the symmetry consid- 
erations above, the current vanishes within the rotating- 
wave approximation. However, as discussed above, we 
expect beyond RWA for small coupling a current / oc F^ . 
Figure |3| confirms this prediction. Yet one observes that 
already a small deviation from </> = is sufficient to re- 
store the usual weak coupling behavior, namely a current 
which is proportional to the coupling strength F. 

The average current for such a harmonic mixing situ- 
ation is depicted in Fig. ^ For large driving amplitudes, 
it is essentially independent of the wire length and, thus. 
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FIG. 3; Average current response to the harmonic mixing 
signal with amplitudes Ai = 2A2 = A, as a function of the 
coupling strength for different phase shifts (j>- The remaining 
parameters are fl = lOA/h, Eb = 5A, ksT = 0.25A. The 
dotted line is proportional to F; it represents a current which 
is proportional to F^. 
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a wire that consists of only a few orbitals, mimics the be- 
havior of an infinite tight-binding system. Figure]^ shows 
the length dependence of the average current for different 
driving strengths. The current saturates as a function of 
the length at a non-zero value. The convergence depends 
on the driving amplitude and is typically reached once 
the number of sites exceeds a value of w 10. For low 
driving amplitudes the current response is more sensitive 
to the wire length. 



C. Rectification in ratciiet-like structures 

A second possibility to realize a finite DC current is to 
preserve the symmetries of the time-dependent part of 



the Hamiltonian by employing a driving field of the form 

a{t) ^ Asm{nt), (48) 



while making the level structure of t 
metric. An example is shown in Fig. pi 



Solecule asym- 
In this molec- 



ular wire model, the inner wire states are arranged in Ng 
groups of three, i.e. N ~ 2 — 3Ng. The levels in each 
such group are shifted by ±Es/2, forming an asymmet- 
ric saw-tooth like structure. 

Figure |^ shows for this model the stationary time- 
averaged current / as a function of the driving ampli- 
tude A. In the limit of a very weak laser field, we find 
/ oc A'^Es, as can be seen from Fig. |[ This behavior is 
expected from symmetry considerations: On one hand. 
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FIG. 6: Level structure of the wire ratchet with N 
atomic sites, i.e., Ng = 2 asymmetric molecular groups. The 
bridge levels are Eb above the donor and acceptor levels and 
are shifted by ±Es/2. 
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FIG. 7; Time-averaged current through a molecular wire 
that consists of Ng bridge units as a function of the driving 
strength A. The bridge parameters are Eb = lOA, Es = A, 
the driving frequency is = 3A/?!,, the coupling to the leads 
is chosen as Vl ~ ~ O.IA/R, and the temperature is 
fcsT — 0.25A. The arrows indicate the driving amplitudes 
used in Fig. ^. 



the asymptotic current must be independent of any ini- 
tial phase of the driving field and therefore is an even 
function of the field amplitude A. On the other hand, / 
vanishes for zero step size Es since then both parity sym- 
metries are restored. The A^-dependence indicates that 
the ratchet effect can only be obtained from a treatment 
beyond linear response. For strong laser fields, we find 
that / is almost independent of the wire length. If the 
driving is moderately strong, / depends in a short wire 
sensitively on the driving amplitude A and the number 
of asymmetric molecular groups Ng] even the sign of the 
current may change with TVg, i.e. we find a current rever- 
sal as a function of the wire length. For long wires that 
comprise five or more wire units (corresponding to 17 or 
more sites), the average current becomes again length- 
independent, as can be observed in Fig. ||. This identifies 
the current reversal as a finite size effect. 

Figure Ru depicts the average current vs. the driving 
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demonstrating the proportionality to A^Es for small driving 
amplitudes. All other parameters are as in Fig. ^. At the dips 
on the right-hand side, the current / changes its sign. 
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FIG. 9: Time-averaged current as a function of the number 
of bridge units Nn, N — 3Ng + 2, for the laser amplitudes 
indicated in Fig. ^. All other parameters are as in Fig. ^. 
The connecting lines serve as a guide to the eye. 



frequency fi, exhibiting resonance peaks as a striking fea- 
ture. Comparison with the quasienergy spectrum reveals 
that each peak corresponds to a non-linear resonance be- 
tween the donor/acceptor and a bridge orbital. While 
the broader peaks at hH. « Eb — lOA match the 1:1 res- 
onance (i.e. the driving frequency equals the energy dif- 
ference), one can identify the sharp peaks for h^l < 7 A 
as multi-photon transitions. Owing to the broken spa- 
tial symmetry of the wire, one expects an asymmetric 
current-voltage characteristic. This is indeed the case as 
depicted with the inset of Fig. ITo|. 
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FIG. 10: Time-averaged current as a function of the angular 
driving frequency Q, for Ng — 1. All other parameters are as 
in Fig. ^. The inset displays the dependence of the average 
current on an externally applied static voltage V, which we 
assume here to drop solely along the molecule. The driving 
frequency and amplitude are f2 = 3A/h (cf. arrow in main 
panel) and A — A, respectively. 



IV. CONCLUSIONS 

With tbis work we have detailed our recently presented 
approachL3 for the computation of the current through 
a time-dependent nanostructure. The Floquet solutions 
of the isolated wire provide a well-adapted basis set that 
keeps the numerical effort for the solution of the master 
equation relatively low. This allows an efficient theo- 
retical treatment that is feasible even for long wires in 
combination with strong laser fields. 

With this formalism we have investigated the possi- 
bility to rectify with a molecular wire an oscillating ex- 
ternal force brought about by laser radiation, thereby 
inducing a non-vanishing average current without any 
net bias. A general requirement for this effect is the 
absence of any reflection symmetry, even in a general- 
ized sense. A most significant difference between "true" 
ratchets and molecular wires studied here is that the lat- 
ter lack the strict spatial periodicity owing to their finite 
length. However, as demonstrated above, already rela- 
tively short wires that consist of approximately 5 to 10 
units can mimic the behavior of an infinite ratchet. If the 
wire is even shorter, we find under certain conditions a 
current reversal as a function of the wire length, i.e. even 
the sign of the current may change. This demonstrates 
that the physics of a coherent quantum ratchet is richer 
than the one of its units, i.e. the combination of coher- 
ently coupled wire units, the driving, and the dissipation 
resulting from the coupling to leads bears new intriguing 
effects. A quantitative analysis of a tight-binding model 
has demonstrated that the resulting currents lie in the 
range of 10~^ Ampere and, thus, can be measured with 
today's techniques. 



An alternative experimental realization of the pre- 
sented results is possible in semiconductor heterostruc- 
tures, where, iastead of a molecule, coherently coupled 
quantum dotsc3 form the central system. A suitable ra- 
diation source that matches the frequency scales in this 
case must operate in the microwave spectral range. 
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APPENDIX A: PARITY OF A SYSTEM UNDER 
DRIVING BY A DIPOLE FIELD 

Although we describe in this work the molecule within 
a tight-binding approximation, it is more convenient to 
study its symmetries as a function of a continuous posi- 
tion and to regard the discrete sites as a special case. Let 
us first consider a Hamiltonian that is an even function 
of X and, thus, is invariant under the parity transforma- 
tion V : X ^ —X. Then, its eigenfunctions ipa can be 
divided into two classes: even and odd ones, according 
to the sign in (fa{x) = ±Lpa{—x). 

Adding a periodically time-dependent dipole force 
xa{t) to such a Hamiltonian evidently breaks parity sym- 
metry since V changes the sign of the interaction with 
the radiation. In a Floquet description, however, we deal 
with states that are functions of both position and time — 
we work in the extended space {x, t\. Instead of the sta- 
tionary Schrodinger equation, we address the eigenvalue 
problem 

H(a;,i)$(a;,t) = e$(x,t) (Al) 
with the so-called Floquet Hamiltonian given by 

H(t) = Ho(x] + xa(t) - ih— , (A2) 

ot 

where we assume a symmetric static part, Ha{x) = 
Hq{—x). Our aim is now to generalize the notion of 
parity to the extended space {x,t} such that the overall 
transformation leaves the Floquet equation ( |Al[ ) invari- 
ant. This can be achieved if the shape of the driving a{t) 
is such that an additional time transformation "repairs" 
the acquired minus sign. We consider two types of trans- 
formation: generalized parity and time-reversal parity. 
Both occur for purely harmonic driving, a{t) — siii{Qt). 
In the following two sections we derive their consequences 
for the Fourier coefficients 

<^k{x)^^ rdte''''''^x,t) (A3) 

of a Floquet states <I>(x,t). 
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1. Generalized parity 



It has been notecE^'El'Lj that a Floquet Hamiltonian of 
the form ( |A2| ) with a(t) = sin(rit) may possess degener- 
ate quasienergies due to its symmetry under the so-called 
generalized parity transformation 
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{X,t) {-X,t + TT/n), 



(A4) 



which consists of spatial parity plus a time shift by 
half a driving period. Thi s symmetry is present in the 
Floquet Hamiltonian (A2), if the driving field obeys 
a{t) = —a{t + ir/fl), since then Scp leaves the Floquet 
equation invariant. Owing to Sqp — 1, we find that 
the corresponding Floquet states are either even or odd, 
SGP^{x,t) = ^{-x,t + Tr/p.) = ±^{x,t). Consequently, 
the Fourier coefhcients (^3) obey the relation 



(A5) 



This transformation represents a further generalization 
of the parity V; we will refer to it as time-inversion par- 
ity since in the literature the term generalized parity is 
mostly used in the context of the transformation ( A4 ) . 



Let us now assume that that the Floq uet Hamiltonian 
is invariant under the transformation (A6), TL{x,t) — 
H*{—x, —t), and that $(a;, t) is a Floquet state, i.e., a so- 
lution of the eigenvalue equation (Al) with quasienergy 
e. Then, ^* {—x, —t) is also a Floquet state with the 
same quasienergy. In the absence of any degeneracy, 
both Floquet states must be identical and, thus, we find 
as a consequence of the time-inversion parity iStp that 
^{x,t) — $*(— X, — t). This is not nessecarily the case 
in the presence of degeneracies, but then we are able to 
choose linear combinations of the (degenerate) Floquet 
states which fulfill the same symmetry relation^ Again 
we are interested in the Fourier decomposition (A3) and 
obtain 



2. Time-inversion parity 



(A7) 



A further symmetry is found if a is an odd function 
of time, a{t) = — a(— t). Then, time inversion transforms 



the Floquet Hamiltonian ( A2) into its complex conjugate 
so that the corresponding symmetry is given by the anti- 
linear transformation 



Stp: ($,x,t) ^ ($*,-x,-i). 



(A6) 



The time-inversion discussed here can be generalized 
by an additional time-shift to read t ^ t^ ~ t. Then, we 
find by the same line of argumentation that ^^(a;) and 
$^(— cc) differ at most by a phase factor. However, for 
convenience one may choose already from the start the 
origin of the time axis such that t^ = 0. 
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